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Recent PHOBOS measurements of the excitation function for the pseudo-rapidity dependence 
of elliptic flow in Au+Au collisions at RHIC, have posed a significant theoretical challenge. Here 
we show that these differential measurements, as well as the RHIC measurements on transverse 
momentum satisfy a universal scaling relation predicted by the Buda-Lund model, based on exact 
solutions of perfect fluid hydrodynamics. We also show that recently found transverse kinetic energy 
scaling of the elliptic flow is a special case of this universal scaling. 


A. Introduction 

One of the unexpected results from experiments at the Relativistic Heavy Ion Collider (RHIC) is the relatively strong 
second harmonic moment of the transverse momentum distribution, referred to as the elliptic flow. Measurements 
of the elliptic flow by the PHENIX, PHOBOS and STAR collaborations (see refs. [H, [2, 3 0, jb. Q) reveal rich 
details in terms of its dependence on particle type, transverse ( p t ) and longitudinal momentum ( 77 ) variables, and on 
the centrality and the bombarding energy of the collision. In the soft transverse momentum region (pt < 2 GeV/c) 
measurements at mid-rapidity are found to be well described by hydrodynamical models [7,0- By contrast, differential 
measurement of the pseudo-rapidity dependence of elliptic flow and its excitation function have resisted several 
attempts at a description in terms of hydrodynamical models (but see their description by the SPHERIO model @,113 
or the approximate descriptions in refs. [llLll2j|). Here we show that these data are consistent with theoretical, analytic 
predictions that are based on perfect fluid hydrodynamics: Fig. |T] demonstrates that the investigated PHOBOS, 
PHENIX and STAR data [I, H,ll, 3 follow the theoretically predicted scaling law. 


B. Perfect fluid hydro picture 


Perfect fluid hydrodynamics is based on local conservation of entropy tr and four-momentum tensor T vtl , 

^(cru^) = 0, 

a ,t» v = 0 , 


(1) 

( 2 ) 

where u M stands for the four-velocity of the matter. The fluid is perfect if the four-momentum tensor is diagonal in 
the local rest frame, 


T» v = {e + p)u> t u v -pcf v . 


( 3 ) 


Here e stands for the local energy density and p for the pressure. These equations are closed by the equation of state, 
which gives the relationship between e, p and a, typically e = up is assumed, where k is either a constant [13 or 
an arbitrary temperature dependent function [141 ] that uses a non-relativistic approximation. Note also, that a bag 
constant can also be introduced, and the e — B = n(p + B) equation of state can be used [Tbl. [Tfil|. 

We focus here on the analytic approach in exploring the consequences of the presence of such perfect fluids in 
high energy heavy ion experiments in Au+Au collisions at RHIC. Such a nonrelativistic exact analytic solution was 
published in ref. [13, while relativistic solutions were published in refs. [H, Gil O • For a detailed discussion on new 
exact relativistic solutions, see ref. fl6[ - 

A tool, that is based on the above listed exact, dynamical hydro solutions, is the Buda-Lund hydro model of refs, 
mill- This hydro model is successful in describing experimental data on single particle spectra and two-particle 
correlations [2fl, [21| ■ The model is defined with the help of its emission function; to take into account the effects of 
long-lived resonances, it utilizes the core-halo model [ 22 ] . 

The elliptic flow is an experimentally measurable observable and is defined as the azimuthal anisotropy or second 
fourier-coefficient of the one-particle momentum distribution N\(j>). The definition of the flow coefficients is: 


v„ = 


/ 0 27r Ai(p) cos [rup)dip 
fo* N i(p)dy> 


( 4 ) 
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where ip is the azimuthal angle of the momentum. This formula returns the elliptic flow v 2 for n = 2. 


C. Universal scaling of the elliptic flow in the Buda-Lund model 


The result for the elliptic flow, that comes directly from a perfect hydro solution is the following simple scaling 
law 00 


v 2 


hiw) 

Io(w) ’ 


(5) 


where I n (z) stands for the modified Bessel function of the second kind, I n (z) = (1/7r) f* exp(ncos(0)) cos(n6)d9. 

Note that this prediction was derived first in 2001 in a non-relativistic perfect hyclrodynamical solution, see eq. (25) 
of ref. 0- In 2003, it has been extended to the relativistic kinematic domain in ref. 0. Ref. [lS| considered 
a relativistic parameterization that included a multitude of known relativistic solutions, such as the Hwa-Bjorken 
solution [H, 0 or other accelerating and Hubble-type of solutions Eli, and interpolated among these. 

The subject of our current investigation is the testing of this scaling law against recent experimental data, but let 
us first discuss this scaling law. 

In the Buda-Lund hydro model, elliptic flow depends only on momentum space anisotropy, but it does not depend 
on the coordinate space anisotropy. This feature of the Buda-Lund model [251 1 is different from azimuthally sensitive 
Blast-wave models [261 ], where v 2 depends also on the coordinate space distribution. 

In section mi we explain the universal scaling variable w in eq. m and show, how w can be determined from 
measurements. In section El we shall subject this relationship to an experimental test also. 


D. Variable of the universal scaling law of v 2 


Looking at eq. m, one sees that the Buda-Lund hydro model predicts [0 a universal scaling: every v 2 measurement 
is predicted to fall on the same scaling curve I\/Iq when plotted against the scaling variable w. This means, that 
v 2 depends on any physical parameter (transverse or longitudinal momentum, mass, center of mass energy, collision 
centrality, type of the colliding nucleus etc.) only through the scaling variable w. This scaling variable is defined by: 


E k 



( 6 ) 


Here Ek is a relativistic generalization of the transverse kinetic energy, defined as 


E k = 


Pt 
2 m t 


with 


(7) 


m t = m t cosh , (8) 

y being the rapidity, Ay the longitudinal expansion of the source, To the central temperature at the freeze-out 
and m* = \/p\ + m 2 the transverse mass. We note, that at mid-rapidity and for a leading order approximation, 
Ek ~ mt — m, which also explains recent development on scaling properties of v 2 by the PHENIX experiment 
at midrapidity 0, 0|. We furthermore note, that parameter Ar/ has recently been dynamically related 0| to 
the acceleration parameter of new exact solutions of relativistic hydrodynamics, where the accelerationless limit 
corresponds to a Bjorken type, flat rapidity distribution and the Ay —> 00 limit. 

The scaling variable w also depends on the parameter T*, which is the effective, rapidity and transverse mass 
dependent slope of the azimuthally averaged single particle spectra, and on the final momentum space eccentricity 
parameter, e. These can be defined [T~tl . fl3 | by the transverse mass and rapidity dependent slope parameters of the 
single particle spectra in the impact parameter (subscript x ) and out of the reaction plane (subscript y ) directions, 
T x and T y , 

1 

Tl 

e 


Ty 


T X - Ty 
T X + Ty 


(9) 


( 10 ) 
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Figure 1: Elliptic flow data of previous plots versus variable w is shown: Data points show the predicted [l!j] universal scaling. 
Small scaling violations at large w values correspond to V 2 (pt) data for pt > 2 GeV. Note, that the error of w was not plotted 
on this plot, but it can be determined from data analysis, and it is on the order of 5-20%. 


which are thus observable quantities. Note also, that e can also be interpreted as a measure of integrated V 2 , and thus 
setting the absolute scale of v- 2 - 

In the Buda-Lund hydro model [3,13 , the rapidity and the transverse mass dependence of the slope parameters 
is given as 


T x 

Ty 


T 0 + rr 
To 


To + m t Y 2 


T 0 + nr 


( 11 ) 

( 12 ) 


Here a 2 = measures the transverse temperature inhomogeneity of the particle emitting source in the transverse 
direction at the mean freeze-out time. 

We note, that each of the kinetic energy term, the effective temperature T* and the eccentricity e are transverse mass 
and rapidity dependent factors. However, for rn t a 2 Tq,T x and T y , hence e and T* become independent of transverse 
mass and rapidity. This saturation of the slope parameters happens only if the temperature is inhomogeneous, ie 
a 2 > 0 . 

The above structure of w, the variable of the universal scaling function of elliptic flow suggests that the transverse 
momentum, rapidity, particle type, centrality, colliding energy, and colliding system dependence of the elliptic flow 
is only apparent in perfect fluid hydrodynamics: a data collapsing behavior sets in and a universal scaling curve 
emerges, which coincides with the ratio of the first and zeroth order modified Bessel functions 00, when V 2 is 
plotted against the scaling variable w. 
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Interesting is furthermore, that the Buda-Lund hydro model also predicts the following universal scaling laws and 
relationships for higher order flows [li|: V 2 n = I n (w) / Iq{ w ) and v 2 n+i — 0. This is to be tested in a later, more 
detailed analysis. 


E. Comparison to experimental data 


We emphasize first, that the scaling variable w is expressed in eq. © in terms of factors, that are in principle 
measurable (however, these factors are not yet determined directly from experimental data). The elliptic flow V 2 is 
also directly measurable. Hence the universal scaling prediction, eq. © can in principle be subjected to a direct 
experimental test. Given the fact that such measurements were not yet published in the literature, we perform an 
indirect testing of the prediction, by determining the relevant parameters of the scaling variable w from an analysis 
of the transverse momentum and rapidity dependence of the elliptic flow in Au+Au collisions at RHIC. 

Transverse momentum dependent elliptic flow data at mid-rapidity can be compared to the Buda-Lund universal 
scaling prediction of 2001 and 2003 of the Buda-Lund model directly, as it was done in e.g. ref. [l9| . 

Eq. (O depends, for a given centrality class, on rapidity y and transverse mass mt . When comparing our result to 
772 ( 77 ) data of the PHOBOS Collaboration, we have performed a saddle point integration in the transverse momentum 
variable and performed a change of variables to the pseudo-rapidity 77 = 0.5 log(j^p^), similarly to ref. f^. This way, 
we have evaluated the single-particle invariant spectra in terms of the variables 77 and <j>, and calculated 772 ( 77 ) from 
this distribution, a procedure corresponding to the PHOBOS measurement described in ref. jlj. 

Scaling implies data collapsing behavior, and also is reflected in a difficulty in extracting the precise values of 
these parameters from elliptic flow measurements: due to the data collapsing behavior, some combinations of these 
fit parameters become relevant, other combinations become irrelevant quantities, that cannot be determined from 
measurements. This is illustrated in Fig. 1, where we compare the universal scaling law of eq. (JSJ) with elliptic flow 
measurements at RHIC. This figure shows an excellent agreement between data and prediction. We may note the 
small scaling violations at largest w values, that correspond to elliptic flow data taken in the transverse momentum 
region of p t > 2 GeV. 

The observed scaling itself shows, that only a few relevant combinations of To, a 2 , X 2 , Y 2 determine the transverse 
momentum dependence of the V 2 measurements. Hence from these measurements it is not possible to reconstruct all 
these four source parameters uniquely. We have chosen the following to eq. © approximative formulas to describe 
the scaling of the elliptic flow: 


2 a 

= jse wr and 

w(pt) = A(1 + B’(m t ~ m) + C’(m t - m) 2 ), 

4 rrit 

and for small values of w eq. © simplifies to V 2 ~ 777 / 2 . The coefficients are as follows: 


A 

B 

A' 


B' 


C' 


Ek 


m t = (m t ),y =0 

1 




m t = (m t ),y =0 


2e 

t! 


mt=m,y —0 

1 T 0 
m Tq T ma 2 


Tn 


1 - 2 - 


T* 


1 


m t —m,y—0 
X 

mt—m,y =0 


l 

m,\T 0 + ma 2 J 
(X 2 + a 2 +Y 2 ) (T 0 + ma 2 ) 3 + 
mX 2 Y 2 (ni 2 (X 2 Y 2 + a 2 (X 2 + F 2 ))) 


- 3 mX 2 Y 2 T 0 (T 0 + ma 2 ) 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 
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V 2 (p) 

20GeV 

62GeV 

130GeV 

200GeV 

A 

0.035±0.004 

0.043±0.001 

0.046T0.001 

0.048±0.001 

B 

0.53T0.1 

0.41T0.01 

0.34T0.01 

0.33T0.01 

X 2 /Ndf 

1.7/11 

9.3/13 

17/15 

18/15 

CL 

91% 

74% 

30% 

28% 

Mv) 

3-15% 

15-25% 

25-50% 


A 

0.028±0.002 

0.048T0.002 

0.061T0.002 

B 

0.64T0.08 

0.60T0.06 

0.43T0.04 

X 2 /Ndf 

12/13 

8/13 

4/13 

CL 

51% 

84% 

96% 


V2 {Pt) 

7T 

K 

P 

A’ [10 4 /MeV] 

5.4T0.1 

6.4T0.3 

3.0T0.1 

B’ [10 4 /MeV] 

16±1 

-0.2T0.4 

25, fixed 

C’ [10 6 /MeV 2 ] 

-1.5T0.1 



X 2 /Ndf 

96/27 

17/5 

27/26 

CL 

lxl0 7 % 

0.5% 

40% 

V2 (pt) 

7T 

K 

P 

A’ [10~ 4 /MeV] 

7.8T0.2 

7.4T0.3 

5.8T0.1 

B’ [10 4 /MeV] 

1.4T0.6 

-1.3T0.4 

1.6, fixed 

C’ [10 7 /MeV 2 ] 

-1.6T0.3 



X 2 /Ndf 

21/10 

13/9 

17/7 

CL 

2% 

15% 

2% 


Table I: Values of the parameters and the quality of the fits for collision energy dependent PHOBOS V 2 (p) data [If is shown in 
the top table, the same for centrality dependent PHOBOS i >2 (p) data [2-f in the second table. The third shows STAR [4jj, the 
fourth PHENIX @] « 2 (pt) data results. 


From this simple picture we had to deviate a little bit in case of proton V 2 (pt) data, here only one parameter could 
have been used to find a valid Minuit [30f minimum, so we fixed B’ there. 

For the case of kaons and protons, only A' and B' were significant, while pion data were so detailed, that C' could 
have been determined, too. Thus we used it only when fitting pion V 2 (pt) • 

For the analysis of the PHOBOS V 2 (p) measurements at RHIC, we have excluded points with large rapidity from 
lower center of mass energies ^(p) fits (77 > 4 for 19.6GeV, p > 4.5 for 62.4GeV). Points with large transverse 
momentum ( p t > 2.0 GeV) were excluded from PHENIX and STAR i >2 (pt) fits. These values give a hint at the 
boundaries of the validity of the model. 

Fits to PHOBOS Ql, PHENIX 1 and STAR Q data are shown in Figs. [2] and [3] The values of the parameters 
and the quality of the fits are summarized in Table HI 

Using the fit parameters in Tableland eqs. (113H14I) . we have determined the universal scaling variable w from these 
PHENIX, PHOBOS and STAR data. 

Using these w values we have plotted the data against w in Fig. [TJ Note that on this plot w itself has an error, as it 
is a reconstructed variable based on eqs. (I15H19D and the values and the errors of the parameters as given in Tabled] 
Standard error propagation has been applied and we have obtained that the relative error of w changes between 
5-50%, and is above 20% only for V 2 (rj) points with large 77 and pion ^(p*) points with large p t . As Fig. |T] is best 
seen when the values of w are plotted on a logarithmic scale, these relative errors of w do not change the qualitative 
conclusion. Comparing with the solid black line in Fig. |T] we observe, that elliptic flow data from various STAR, 
PHENIX and PHOBOS measurements at RHIC follow the universal scaling curve, predicted by the ellipsoidally 
symmetric Buda-Lund model in 2004. 


I. FURTHER SCALING PROPERTIES 

From the Taylor expansion of the Bessel functions for the realistic w <C 1 case one finds 

w 



(20) 
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Figure 2: PHENIX Q and STAR [j data on elliptic flow, V 2 , plotted versus pt and fitted with Buda-Lund model. 


One can also see, that for small momenta, 


E k = 7 —^— « m t - m. ( 21 ) 

2 m t 

Thus a leading order calculation, at mid-rapidity, from eq. m one gets 

A' 

V2 «—( 22 ) 

This derivation indicates, that the PHENIX discovery [27j, [28| of the universal scaling of the elliptic flow in terms of 
mt — m at mid-rapidity is a consequence, a special case of the more general universal scaling law of eq. [5j predicted 
by the Buda-Lund model. 

Furthermore, PHENIX found [ 27 ], [28|, that V 2 /n q scales in terms of (— m)/n q for several types of particles ( 7 r, K, 
K°, p, A, S, deuteron, d>). This indicates that the perfect fluid motion scales on the quark level, and the Buda-Lund 
model scaling prediction for higher order flows, V 2 n = I n {w)/Io(w) (see ref. 0 ) should also be applied on the quark 
level. 









n 





Figure 3: PHOBOS 0,1 data on elliptic flow, V 2 , plotted versus rj and fitted with Buda-Lund model. 


A. Conclusions 

We have shown that the excitation function of the transverse momentum and pseudorapidity dependence of the 
elliptic flow in Au+Au collisions is well described with the formulas that are predicted by the Buda-Lund type of 
hydrodynamical calculations 00. We have provided a positive test for the validity of the perfect fluid picture of 
soft particle production in Au+Au collisions at RHIC up to 1-1.5 GeV and up to a pseudorapidity of rybeam — 0.5. 

We have also shown that the PHENIX discovery [27), (28| of a scaling behavior of v% vs. mt — m is a special case 
of the more general, rapidity dependent universal scaling law of the Buda-Lund type of perfect fluid hydrodynamical 
solutions. 

The universal scaling of PHOBOS V 2 (tj) and PHENIX and STAR v 2 (pt) , expressed by eq. © and illustrated by 
Fig. []] provides a successful quantitative as well as qualitative test for the appearance of a perfect fluid in Au+Au 
collisions at various colliding energies at RHIC. 

We have furthermore shown, that since PHENIX found 03, 0, that V 2 /n q scales in terms of (mt — m)/n q for 
several types of particles, the Buda-Lund model scaling prediction for higher order flows, V 2 n = In(w) / Iq(w) (see 
ref. pj|) should also be applied on the quark level. 
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